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ABSTRACT 

We have applied Pade approximants to perturbative QCD calculations of event shape 
observables in e + e~ — > hadrons. We used the exact O(a^) prediction and the [0/1] 
Pade approximant to estimate the O(ctg) term for 15 observables, and in each case 
determined a s (M^) from comparison with hadronic Z® decay data from the SLD 
experiment. We found the scatter among the a s (M^) values to be significantly reduced 
compared with the standard O(ttg) determination, implying that the Pade method 
provides at least a partial approximation of higher-order perturbative contributions to 
event shape observables. 
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One of the most important tasks in high energy physics is the precise 
determination of the strong coupling a s , conventionally expressed at the scale of the 
mass of the Z° boson, Mz — 91.2 GeV. Not only does measurement of ck s (M|) in 
different hard processes and at different hard scales Q provide a fundamental test 
of the theory of strong interactions, Quantum Chromodynamics (QCD), but it also 
allows constraints on extensions to the Standard Model of elementary particles. For 
example, it has been claimed [1] that measurements of electroweak observables can 
be better described by the Standard Model with a s (Mg) = 0.112 and the addition 
of light superpartners, than by the Standard Model alone with a s (M|) = 0.123. The 
latter study, as well as studies of other possible extensions to the Standard Model 
[2], have been prompted by claims [3] that recent measurements of a s (M|) may be 
grouped into two classes: those made at 'low-Q 2 ', which tend to cluster at values 
around 0.112, and those made at 'high-Q 2 ', which tend to cluster at values around 
0.123; for a review of these measurements see Ref. [4]. 

Examination of the large set of a s (M|) measurements reveals that in fact all are 
consistent with a 'world average' central value of about 0.117 with an uncertainty 
of ±0.005 [4], and that their grouping into two supposedly discrepant classes is 
arbitrary and not significant. Furthermore, nearly all measurements are limited by 
theoretical systematic uncertainties that derive from lack of knowledge of higher-order 
perturbative QCD contributions, or of non-perturbative effects, or both. The effects 
on a s (M|) determinations of such uncalculated contributions have been estimated 
using ad hoc procedures, sometimes in different ways by different experimental 
collaborations working in similar areas; see eg. [5,6]. It appears to us premature 
to speculate on beyond-Standard Model explanations of a s (M|) measurements 
that, within large, dominant and arbitrarily-estimated theoretical uncertainties, are 
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consistent with one another. We suggest that a more rational approach is to reduce 
the size of the limiting theoretical uncertainties which may, or may not, be concealing 
new physics. 

Here we consider hadronic event shape observables in e + e~ annihilation, for 
which perturbative QCD calculations have been performed exactly only up to 
second order in a s [7,8] and have been used extensively by collaborations at the 
PETRA, PEP, TRISTAN, SLC, and LEP colliders for measurement of a s (M 2 z ) [4]. 
A large scatter among the a s (M z ) values derived from different observables has 
been found [9,10], which can in principle be accounted for as an effect of the a 
priori unknown higher order contributions in perturbation theory. A consensus 
has arisen among experimentalists to estimate the size of this effect for each 
observable from the renormalisation-scale (fx) dependence of the a s (M z ) values 
derived from fits of the 0(a 2 ) calculation to the data, see eg. [10], and to 
quote a corresponding renormalisation scale uncertainty on a s (Ml). Within such 
uncertainties the a s (M z ) values from the different observables are found to be 
consistent, but this (arbitrary) procedure naturally leads to a large uncertainty on 
the final a s (M z ) averaged over all observables; for example, using 15 observables the 
SLD Collaboration determined [10]: 

a s (M\) = 0.1226 ± 0.0025 (exp.) ± 0.0109 (theor.) (1) 

where the theoretical uncertainty is dominated by the contribution of ±0.0106 from 
the renormalisation scale variation. 

The best resolution of this situation would be to calculate the observables to 
higher order in perturbation theory, a difficult and unattractive task that has not 
yet been achieved. In the absence of 0(0;^) QCD calculations it has been suggested 
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[11,12,13] that the O(a^) calculation for each observable can be 'optimised' by 
choosing a specific value of the renormalisation scale. It has recently been shown 
[6] that such optimised scale choices do not serve to reduce the scatter among the 
a s (Mz) values determined from different observables, which implies that they do not 
reduce the size of the contributions from the uncalculated higher orders. 

Here we employ an approach for estimating the 0(otg) contributions to, as well 
as the sum of, the perturbative QCD series for event shape observables in e + e~ 
annihilation. The method makes use of Pade Approximants (PA). The PA method 
has been applied to estimate coefficients in perturbative quantum field theory and 
statistical physics [14] and is outlined in detail in Ref. [15]. We give a brief review 
here. 

The PA [N/M] to the series: 

S = S + S lX + S 2 x 2 + ... + S N+M x N+M (2) 

is defined [16]: 

= a + aix + a 2 x z + ... + a N x J " 
[ 1 J ~ 1 + hx + b 2 x 2 + ... + b M x M ' 1 ' 

where N and M are integers such that N > and M > 0, and 

[N/M] = S + 0(x N+M+1 ). (4) 

The coefficients a{ (0 < i < N) and bj (1 < j < M) are obtained by multiplying 
Eq. 4 by the denominator of Eq. 3 and equating coefficients of like powers of x. By 
consideration of the terms of 0(x N+M+1 ) one can obtain an estimate of the coefficient 
Sn+m+i- Furthermore, for an asymptotic series [N/M] can be taken to be an estimate 
of the sum of the series to all orders; we refer to this as the Pade sum (PS). 
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Recently this method has been applied to estimate the O(af) term in 
the perturbative QCD prediction for the inclusive quantities R = a(e + e~ — > 
hadrons)/o-(e + e~ — > /i + p~), R T = T(t^ v hadrons)/ T(r^ eV e v T ) and the Bjorken 
sum rule in deep inelastic scattering [17]. In each case the PA estimate is in good 
agreement with the exact calculation of the O(of) term, which is remarkable since 
the method does not involve knowledge of strong interaction dynamics. Based upon 
this success PA estimates have also been made of O(af) terms for the same quantities 
[17,18]. 

In the case of e + e~ annihilation into hadronic final states the perturbative QCD 
prediction for an infra-red- and collinear-safe observable X can be written: 

-%{X,n) = a- s (fi)A(X) + af{^) B{X^) + af{^) C(X, fi) + O(of(fi)) (5) 

where = a s /2n. To date only the leading and next-to-leading coefficients A(X) and 
B(X, n) have been calculated [7,8]. Throughout this paper we set the renormalisation 
scale fj, to the so-called physical value /i — Q — y/s; in this case the //-dependence of 
the beyond-leading coefficients B, C, . . . vanishes. 

We have employed the 15 hadronic event shape observables used in the recent 
a s (Mz) determination by the SLD Collaboration [10], namely 1— thrust (r), 
oblateness (O), the C-parameter (C), normalised heavy jet mass (p), total jet 
broadening (Bt), wide jet broadening (Byy), the differential 2-jet rate defined using 
the E, E0, P, P0, D and G algorithms (D%, £>f°, Dg, D^°, D§ and D% respectively), 
energy-energy correlations (EEC) and their asymmetry (AEEC), and the Jet Cone 
Energy Fraction (JCEF) [19]. Distributions of these 15 event shape observables were 
measured [10] using a sample of approximately 50,000 hadronic Z° decay events. 
The data were corrected for detector bias effects such as acceptance, resolution, and 



inefficiency, as well as for the effects of initial-state radiation and hadronisation, to 
arrive at 'parton-level' distributions, which can be compared directly with the QCD 
calculations. 

For each observable X we employed the EVENT program [20] to calculate the 
coefficients A(X) and B(X) in Eq. 5. These coefficients are listed in Table 1 for 
a representative value of each observable. In each case it can be seen that the 
next-to-leading coefficient B is typically an order of magnitude larger than the leading 
coefficient A. At first sight this appears to make the perturbative approach invalid; it 
should be noted, however, that the ratio of next-to-leading to leading terms in Eq. 5 
contains an additional factor of aj, which is about 1/52 for a s (M§) = 0.12, so that 
in all cases the next-to-leading term is smaller than the leading term. 

By explicitly separating an overall factor of on the r.h.s. of Eq. 5 and comparing 
with the form of Eq. 2 the PA [0/1] can be defined for the series. For each bin of each 
observable, by consideration of the PA [0/1] we derived an estimate of the coefficient 
C(X) of the next-to- next-to-leading term in the series: 

dX) = igf (6) 

Examples of the predictions for C(X) are given in Table 1. It follows from Eq. 6 that 
C(X) > and that the ratios C{X)/B(X) and B(X)/A{X) are equal. For each 
bin of each observable we added the PA prediction for the O(al) term to the exact 
O(a^) calculation to obtain an estimate of the series to O(a^). For each observable 
we then fitted these calculations simultaneously to all bins in the selected range of 
the data [10] by minimising x 2 w.r.t. variation of a s (M|). The resulting a s (M|) and 
Xdof vames are shown in Table 2. 
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It should be noted that the 0(a^) estimate does not provide a good fit to the Bt 
data, as indicated by the large Xdof varue (Table 2). It is perhaps significant that this 
observable shows a large renormalisation scale uncertainty on a s (M|) determined 
at O(a^) [10], and has a very large ratio of next-to-leading to leading order 
coefficients (Table 1), both implying that beyond- next-to-leading order contributions 
may be large. Furthermore, the latter suggests that the PA [0/1] estimate of the 
next-to-next-to-leading order coefficient may be neither reliable, nor sufficient to 
describe the data as is observed. In the following discussion we therefore exclude 
Bt and restrict the analysis to the remaining 14 observables; however, we indicate in 
footnotes to the text any relevant changes in results if Bt is included. 

The a s (Mz) values are shown in Fig. 1, together with the corresponding values 
from the O(a^) analysis of the SLD data with \x = y/s [6]. For each observable it 
can be seen that the a s (M|) value derived using the 0(ofy estimate is lower than 
that derived using the O(a^) calculation, which is expected since C(X) is positive. 
Also shown for each observable is the total uncertainty on a s (M|) from the SLD 
0(a!g) study [10], which is dominated by the renormalisation scale variation. In each 
case the O(a^) a s (M^) value lies near the lower bound given by the scale uncertainty 
on the 0(a>g) result. To the extent that the PA 0(af ) estimate is accurate, this implies 
that the uncertainty assigned to the O(a^) a s (M^) value from each observable due 
to variation of the renormalisation scale is a reasonable estimate of the effect of the 
missing 0(dg) contribution. 
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Furthermore, it can be seen from Fig. 1 that the scatter among the a s (M z ) values 
is noticeably smaller in the 0(a^) case. Taking an unweighted average 1 and r.m.s. 
deviation over each set of 14 a s (M z ) values yields : 

a s (M z ) = 0.1255 ± 0.0070 (0(a 2 s )) (7) 

a s {M 2 z ) = 0.1147 ± 0.0035 (0(a 3 s )), (8) 

implying that the inclusion of 0(0^) terms causes a reduction in the central value 
of a s (M z ) at /i = y/s by approximately 0.011, and that the residual scatter 
among the O(af) a s (M z ) values, presumably due to missing O(ag) terms, is 
approximately ±0.0035, which is comparable with the combined experimental error 
and hadronisation uncertainty on a single observable measured at Q = Mz [10]. Since 
the accuracy of the Pade Approximant method can only be verified a posteriori, exact 
calculation of the O(a^) terms in order to confirm these results would be extremely 
desirable. 

We also used the PS [0/1] as an estimate of the sum of the asymptotic series 
and extracted a s (M z ) by comparison with the data in a similar manner. The 
a s (M z ) values are shown in Fig. 1. Typically, for each observable, the PS 
a s (M z ) value is close to the O(a^) value. Again the fit to Bt is very poor, and 
we omit it from the following main discussion. Taking an unweighted average and 
r.m.s. deviation over the set of 14 a s (M z ) values yields : 

a s {M 2 z ) = 0.1148 ± 0.0052 (PS). (9) 

1 Weighted averages using experimental errors [10] yield mean a s (M|) values consistent with 
the respective unweighted averages within the statistical error on a single a s {M\ ) value. 

2 Including Bt and averaging over 15 observables yields a s (M\ ) = 0.1265 ± 0.0076 (0(a 2 )) 
and a s (M|) = 0.1139 ± 0.0045 (0(a 3 s )). 

3 Including Bt and averaging over 15 observables yields a s (M|) = 0.1136 ± 0.0068. 
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Though the average value is close to that obtained using the PA O(a^) estimate, the 
r.m.s. deviation is somewhat larger, implying that the PS [0/1] provides a poorer 
estimate of the sum of the series than the PA [0/1] estimate to O(a^). 

Other interesting features are apparent from Fig. I. Several observables yield 
noticeably larger scale uncertainty at 0(cig), namely r, Bt, O, C, and EEC, 
which is an indication of potentially large contributions at O(a^) that is supported by 
the PA and PS [0/1] estimates. If these are omitted from consideration the remaining 
set of observables yields average and r.m.s. a s (M§) values of 0.1212 ± 0.0044 
(0{a 2 s ) analysis), 0.1131 ± 0.0028 (0(af) analysis) and 0.1155 ± 0.0025 (PS analysis). 
Though the selection of this subset is arbitrary, it is perhaps noteworthy that both 
Pade-derived a s (M|) values are close to the corresponding values [6] obtained with 
0(a 2 s ) calculations and 'PMS' [11] and 'FAC [12] optimised scales, 0.1146 ± 0.0025 

4 

and 0.1148 ± 0.0025 respectively, and that the scatter is similarly small . 

In summary, we have applied Pade approximants to the determination of 
a s (M§) from hadronic event shape observables in e + e _ annihilation. We applied 
the PA [0/1] to the 0(a 2 ) perturbative QCD series for 15 observables to obtain 
estimates of the O(af) coefficients, and then fitted the extended series to SLD data 
to determine a s (M|). With the renormalisation scale fixed to the cm. energy 
the scatter (neglecting Bt) among the a s [M\) values was reduced from ±0.0070 
{0{a 2 s )) to ±0.0035 (0(a 3 s )), and the central value of a s (M§) was lowered by 
0.011. If the scatter is interpreted as arising from missing higher-order contributions, 
the reduction in scatter implies that the Pade method provides at least a partial 
approximation of higher-order perturbative QCD contributions to event shape 
observables. Furthermore, if the Pade-estimated 0(a 3 ) terms are accurate, this result 

4 Significantly larger scatter is obtained with the 'BLM' scale choice: a s (M§) = 0.1082 ±0.0091. 
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implies that residual 0(af) terms contribute to c* s (M|) at the level of ±0.0035, which 
is comparable with current experimental and hadronisation uncertainties. 

These results based on Pade approximants are tantalising, but they can only 
be verified upon completion of a full perturbative QCD calculation at 0(a!g), which 
we strongly encourage. One could then apply the [0/2] or [1/1] Pade approximants 
to the 0(ag) series in order to estimate the size of 0(af) contributions. Since the 
accuracy of the Pade method in predicting unknown perturbation series coefficients 
is expected to improve with increasing order of the known terms, it can be argued 
that the estimated O(af) coefficients would be expected to be more accurate than 
the estimates of O(otg) coefficients presented here 5 . Combined with recent theoretical 
progress in understanding hadronisation effects in terms of 'power corrections' [21], 
it may then be plausible to expect that a s (M§) measurements at the 1%-level of 
precision could be achieved at future high-energy e + e~ colliders [22] . 

We thank D. Muller and Y. Ohnishi for helpful contributions, and also S. Brodsky, 
L. Dixon and T. Rizzo for comments on the manuscript. P.N.B. thanks Lawrence 
Morland for useful discussions concerning Pade approximants. 



5 See Refs. [17,18] for O(a^) estimates for R, R T and the Bjorken sum rule. 
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Observable 


1 ~\7~ 

value X 


A(X) 


B(X) 


B{X)/A{X) 


C(X) 


: t 


0.18 


23.93 


711.1 


29.7 


21130 


P 


0.15 


36.55 


624.5 


17.1 


10670 


Bf 


0.23 


21.69 


1042 


48.0 


50020 


B w 


0.14 


82.16 


1402 


17.1 


23940 





0.21 


71.27 


-363.0 


-5.09 


1848 


c 


0.46 


18.89 


518.9 


27.8 


14260 


Vc (Df) 


0.17 


38.84 


1217 


31.3 


38140 


Vc (Df u ) 


0.17 


38.84 


822.1 


21.2 


17400 


Vc (Df) 


0.12 


78.10 


1029 


13.2 


1 or rn 

13550 


Vc (£>2 ) 


0.17 


OO OA 

38.84 


601.6 


15.5 


r\ "i 

9318 


Vc {D§) 


0.12 


33.10 


490.6 


14.8 


7271 


Vc (D°) 


0.22 


37.93 


412.7 


10.9 


4489 


X (EEC) 


91.8° 


2.460 


44.23 


18.0 


795 


X (AEEC) 


41.4° 


2.682 


22.74 


8.48 


193 


X (JCEF) 


131.4° 


5.196 


50.97 


9.81 


500 



Table 1. Perturbative QCD calculation of 1/ada/dX for each observable (first 
column) at a representative value X (second column). The third and fourth columns 
show the leading and next-to-leading order coefficients, respectively. The ratio 
of coefficients is shown in the fifth column, and the PA [0/1] prediction for the 
next-to-next-to- leading order coefficient C(X) is listed in the sixth column. 
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Observable 


a s {Ml) 


2 

Xdof 


T 





1172±0.0009 


0.6 


P 





1183±0.0007 


0.9 







1031±0.0009 


112 


B w 





1143±0.0008 


1.6 








1230±0.0011 


0.6 


C 





1158±0.0008 


2.7 


D 2 





1164±0.0006 


6.4 


U 2 





1105±0.0007 


3.5 


7-1 P 





1144±0.0008 


0.8 


nPO 





1122±0.0009 


3.4 


u 2 





1156±0.0011 


4.0 


D G 

2 





1122±0.0008 


3.9 


EEC 





1154±0.0008 


0.4 


AEEC 





1082±0.0012 


0.4 


\ JCEF 





1124±0.0007 


0.2 



Table 2. Values of a s {M\) and Xdof fr° m fits of QCD predictions incorporating the 
PA [0/1] estimate for the O(a^) term. The errors are statistical only and are highly 
correlated between observables. 
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Figure Caption 



FIG. 1. Values of a s (M|) determined from fits to event shape observables (see text): 
(a) O(aig) (circles); (b) 0(a^) estimate (squares); (c) Pade sum (PS) (crosses). The 
shaded band shown for each observable is dominated by the renormalisation scale 
uncertainty on the 0(af) fit. For each point combined experimental statistical and 
systematic errors [10] are shown; the errors are highly correlated between observables. 
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